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W) ABSTRACT 

This paper defines the beta function and other linear orbit parameters 
,_i using the exact equations of motion. The /3, a and ip functions are redefined 

CO using the exact equations. Expressions are found for the transfer matrix 

and the emittance. The differential equation for rj = x/fi 1 / 2 is found. New 

relationships between a, (3, ij) and v are derived. 
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c-j 1. Introduction 

O This paper defines the beta function and the other linear orbit parameters using the 

^ exact equations of motion. The usual treatment 1 of the linear orbit parameters is based 



on the approximate equation of motion 



d 2 



x 



+ K(s)x = (1.1) 



ds 2 

Approximations are made in obtaining Eq. (1.1) which are usually valid for large acceler- 
ators. 

The exact linearized equations of motion can be written as 

dx 

— = A n x + Ai2p x 

dt W 
— — = A 2 \x + A 2 2Px 

as 
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x and p x are the canonical coordinates in in a curvilinear coordinate system based on a 
reference orbit and the Aij(s) are periodic in s with period L. The approximate Eq. (1.1) 
assumes that A\\ = A22 = 0, A\2 = 1 and A21 = —K(s). The exact values of the A{j are 
given in section 2. 

A treatment of the linear orbit parameters based on the exact equations, Eqs. (1.2), 
rather than the approximate Eq. (1.1) may be desirable in the following situations: 

1. Symplectic long term tracking using a procedure where the magnets are re- 
placed by a sequence of point magnets and drift spaces. For the tracking to be 
symplectic, one has to use the solutions of the exact equations of motion. The 
linearized equations of motion then have the form of Eq. (1.2). 

2. Small accelerators where the approximation made in deriving Eq. (1.1) may 
not be valid. 

Many of the results found using the approximate equations carry over for the exact 
equations. A few of the changed results are the following: 

1 f ld ( 3 AO 

a =M-2\-2Ts +A ^ 

* = I An j (1.3) 

v = — / ds- 12 



2tt Jo (3 
where C is the circumference of the accelerator. 

Some unchanged results are 

e = 7a; 2 + 2«7p + (3p 2 

cos /i + a sin /i (3 sin /i 



M(s + L,s) 

7 sin n cos n — a sin // j (14) 
7= (l + a 2 )/(3 

li = ip (s + L) - ip (s) 

L is the period of the Aij(s). More detailed results are given below. Because of Eqs. (1.4), 
the usual procedure used in tracking programs to compute u } /3, a, 7 and ip, from the 
transfer matrix is still valid. 
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2. Equations of Motion 



The exact equations of motion can be written 1,2 as 

dx 1 + x I p 

T = 

as q s 



dq x 
ds 
dy_ 
ds 
dty 
ds 



q s i 

P Bp 

1 + x/p 



1 + x/p) B y -^- (1 + x/p) B s 

q s 



{2.W 



— (1 + x/p) B a - (1 + x/p) B x 



Bp [q s 

Qs = (l - ql - q 2 y ) 112 , Bp = pc/e 
In Eq. (2.1), the particle coordinates in a curvilinear coordinate system based 

on a reference orbit with the radius curvature p(s). p X} p S} p y are the components of the 
momentum and 

q x = Px/p, q y = Py/p, q s = p s /p (2.2) 

For large accelerators, q x ~ x' and q y ~ y' . 

To find the linearized equations for the betatron oscillations, one expands Eq. (2.1) 
around the closed orbit for a particular momentum, p. This gives the set of linear equations 

—j— = AijXj (2.3a) 
j 

where ~X{ are the 4 coordinates relative to the closed orbit. The A{j are given by 

Aij = ^-fi i,j = l,4 (2.3b) 

OXj 

evaluated on the closed orbit and Eqs. (2.1) have been written as 

fi i = l,4 (2.3c) 



ds 



For the case where B s = 0, q y = 0, y = on the closed orbit, the A{j are given by 



A, 



1 

B~p 



lq^ 

p q s 

B v dB v 

y + (1 + x/p) 



P 



da 








+ 

ql 

_\qx_ 
pq s 
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1 n -l. / \ dB * 

-B- P {1 + X/P) ^y- 







q 3 s 



(2.3d) 
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For large accelerators where g z < 1, 5j < 1, < 1 one has 

An = A 22 = 0, A n = 1 

(2.4) 

Am = 1. 

For the exact equations, An and A22 are not zero, and one does not have A12 = 1, A34 = 1. 
In particular, 

An = -A22 = (1/p) (q x /q s ) 

Ai 2 = (l + x/p)/q 3 s (2.5) 
A M = {l + x/p)(l-q 2 x ) /q 3 s 



3. Eigenfunctions of the Exact Linear Equations of Motion and the Linear Or- 
bit Parameters 

The problem now is, given the exact linear equations of motion, Eq. (2.3), how does 
one define the linear orbit parameters /3, a, 7, v and the emittance e, and what are the 
relationships that hold between them. To do this, one has to repeat the well known 
treatment of the linear orbit parameters, and see where the definitions and relationships 
change for the exact equations. The treatment given below is believed to reduce the amount 
of algebraic manipulation required, and makes few assumptions about the A{j coefficients 
in the linear equations. 

For the x motion, the linear equations are written as 

dx 

— = Anx + Ai2p x 

—— = A 2 ip x + A22X 

as 

The A{j are given by Eqs. (2.3). 

The transfer matrix M(s } sq) obeys 

x = M (s, so) xq 

x 



Px 1 i3.2) 



^—M = AM 

ds 
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One may note that the symbol x is used in 2 different ways. The meaning of x should be 
clear from the context. The matrix M is symplectic as the equations of motion are derived 
from a hamiltonian. 1,2 Thus 

MM = 1 

M=SMS 

S= (-l o) ' /= (o 1 

S is the transverse of S. Also \M\ = 1; \M\ is the determinant of M. 
The one period transfer matrix is defined by 

M(s) = M(s + L,s) (3.4a) 

where L is the period of the A{ y in Eq. (3.1). One can show that M(s) and M(so) are 
related by 

M(s) = M(s,so)M(so)M(sq,s) (3.4b) 

The eigenfunctions and eigenvalues of M(s) are defined by 

M (s) x = \x } 
\M-\I\ = } (3.5) 
A 2 - (m n + m 22 ) A + 1 = 

where mij are the elements of M, and using \M\ = 1. 

Eqs. (3.5) shows that the two eigenvalues Ai, A 2 obey 

AiA 2 = 1, (3.6a) 

and for stable motion, |A| = 1 and A 2 = A*, and we can write 

Ai = exp (i/i) (3.6b) 

Given the eigenfunction at so, ^i(-so) one can find the eigenfunction at any other point 
s using 

x\ (s) = M (s, s ) xi (s ) , (3.7a) 
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and xi(s) has the same eigenvalue Ai. This follows from Eq. (3.5), using Eq. (3.4b) to 
relate M(s) and M(sq). Also xi(s) obeys the linear equations of motion, 



d_ 

as 

which follows from Eq. (3.7a) and Eq. (3.2). One can show that 

xi (s) /Xl /L = fi (s) 
where f\(s + L) = fi(s). This follows from 

f 1 (s + L) = x 1 (s + L)/\' 1 /L+1 



(3.7b) 



(3.8a) 



M(s) Xl (s) /Xl /L+1 = Xl (s) /Xl /L 



Thus, one can write 



(3.8b) 



(3.9) 



x\ (s) = exp (i/is/L) fi (s) 

h (s + L) = h (s) 
Eq. (3.8b) can be rewritten as 

x l ( s ) = ( s ) 1/ ' 2 ex P {i^) 
if; (s) = jjs/L + gi (s) 

gi (s + L)= gi (s), (3 (s + L) = (3 (s) 
Eq. (3.9) defines the beta functions, (3(s), except for a normalization multiplyer, for the 
eigenfunction xi(s). The normalization multiplyer will be defined below. It will be shown 
first that ip and (3 are related. To find this relation, one uses the Lagrange invariant 1 

W=x 2 Sxi (3.10) 

where x\ } x 2 are two solutions of the equations of motion. Eq. (3.10) corresponds to the 
Wronskian in the treatment of the approximate equations of motion. For x\ and x 2} we 
use the two eigenfunctions x\ and x 2 = x*. 



X\ 



Xl 
Pxl 



For xi one uses Eq. (3.9) and for p x i one finds from the equations of motion 

1 / dxi 



Pxl 



Anx 



A12 v ds 

W = X 2 p x l - Px2Xl 

dxi xidx 2 



W 
W 



x 2 - 



ds 



ds 



A 



(3.11) 



(3.121 



(3.13) 



12 



2i ^ 

A\2 ds 
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The beta function (3 is normalized by normalizing the eigenfunctions so that 

W=x* 1 Sx 1 = 2i (3.14) 



which gives 



dip 
ds 



A 



12 



P 



(3.15) 



Eq. (3.15) replaces the familiar result dip / 'ds = 1/(3 which is obtained when A\2 = 1. From 
Eq. (3.15) one can find a result for the tune. Using 2ttv = ip(C) — ip(0) where C is the 



circumference of the accelerator, one finds 



C 



v = — / ds — 

2tt J (3 



(3.16) 



^From Eq. (3.12) we now find for p x \. 



Pxl = ~ «) ex P(^) 



(3.17a) 
(3.17b) 



Eqs. (3.17) provide the new definition for the a parameter, which replaces the familiar 
result a = ^d(3/ds. At this point the definition of a may seem arbitrary. It will be seen to 
be the convenient definition of a when the emittance and transfer matrix are considered 
below. 

The eigenfunctions can now be written as, using Eq. (3.9) and Eq. (3.17), 



{-a + i). 



exp (iip) 



(3.18) 
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4. The Transfer Matrix, the Emittance and the Linear Orbit Parameters 



The particle motion can be written as a linear combination of the eigenfunction given 
by Eq. (3.18) 

x = ax\ + c.c. 

x = 2\a\(3 1/2 cos(V> + 8) 
ax 2\a\ sin 
Px ~~ (3 ~ /?i/2 
a = \a\ exp (iS) 

Eqs. (4.1) suggest the new variables r/ } p v 



Pv 



G 



Px 



G 



/3~2 



for then 



77 = 2|«| cos (ip + 8) 

Pn = — 2 1 a I sin (ip + 8) 
one obtains the emittance invariant from Eq. (4.3) 

2 1 2 

r] + p v = e 



(4.2) 



(4.3) 



(4.4) 



2\a\ = e 1 ' 2 

r] = e 1 / 2 cos (ip + 8) , p n = — e 1//2 sin (ip + 8) 
Replacing r/ } p^, by x, p x using Eq. (4.2), one gets 

e = -fx 2 + 2axp x + [3p 2 x 
7= (l + « 2 ) I [3 

It will be shown that e is the phase space area, divided by 7r, inside the ellipse defined by 
Eq. (4.5). 

Since |Cr| = 1, G is a symplectic matrix and 77, p v are symplected variables with the 
transfer matrix U(s,sq) and 



(4.5) 



Pv 



U (s,s ) 



Pvo 



(4.6) 



U = G(s)M(s,s )G(s c 
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G = G -1 



fj 1 / 2 



U = U(s + L, s) has the same eigenvalues as M, and the eigenfunction of U, r/i and 
= ^1 are related to the eigenfunctions of M, x\ and £2 = £*, by 



771 = Gxi 
1 



^1 



(4.7) 



where Eq. (3.18) was used for x\. 

One sees from Eqs. (4.3) that p v = drj/dtp and thus Eq. (4.3) can be rewritten as 

?7 = 770 cos (ip - ip ) + P-qO sin (ip - ipo) 
P v = -r]o sin (ip - V'o) + P v cos (V> - V>o) 
Eq. (4.8) gives a result for U(s,sq) 



(4.8) 



U (s, Sr 



cos (tp — ipo) sin [ip — ipo) 
— sin (ip — ipo) cos (?/> — ?/>o) 



One can then find M(s } sq) using M = G(s)UG(sq) } 
1/2 



M{s,s ] 



One then finds 



[cos (ip - ipo) + «o sin (?/> - ipo)] 



'1 + aao 



sin(^-^o) 



Wo) 



1/2 



M = M (s + L,s) ■ 
= ip (s + L) - ip (s) 



cos /j + a sin /i 
—7 sin /i 



1/2 



(4.9) 



(/?/?o) 1/2 sin(V>-Vo) 

cos (ip — ipo) — a sin (?/> — ipo)] 
(4.10) 



/3 sin /i 
cos /i — a sin /i 



(4.11) 



The results for M and M are unchanged from those found for the approximate equations of 
motion. These results and the result for the emittance Eq. (4.5) justify the new definition 
for a, Eq. (3.17). 

The connection between the emittance, defined by Eq. (4.5), and the phase space area 
inside the ellipse defined by Eq. (4.5), is given by 



J dxdp x = J drjdpjj = ne 



(4.12) 



where we have used \G\ = 1 and Eq. (4.4) for the ellipse in 17 , p v space. 
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The relationship between the parameters (3, a, 7 at s and (3, a, 7 at sq is unchanged. 



To see this, write 



M = I cos u + J sin it 

a (3 
—7 —a 



J 



Using Eq. (3.4b) that connects M(s) and M(sq) one finds 

J = M(s,so)J(so)M- 1 (s,s ] 
Eq. (4.13) gives the desired result 











.7. 





m 



— 2m\m2 
-1112111111 1 + 2mi2m2i 
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m 



21 



-2m22"^21 



m 12 

-mi2m22 
,2 



m 



22 













_7o _ 



(4.13) 



(4.14) 



(4.15) 



5. Differential Equations for the Linear Orbit Parameter 

This section finds differential equations for (3, r/ and (3, a, 7. 



5.1. Second Order Differential Equation for x 

^From the first order differential equation for x, p X} Eq. (3.1), one can eliminate p x to 
find a second order equation for x. ^From Eq. (3.1) 





1 


{ dx 


p = 




\ ds 


dp 


d 


1 


ds 


ds 


A 12 


d 


( 1 


ds \ 


ds 


\AndsJ 



- A\\x 
dx 

A\\x 

as 

+ x [ -A21 



Mix + — — ( Anx 



d f Au 



ds\A 



L12 



An \ ds 



A 2 



A 



12 



(5.1) 



It has been assumed that A\\ = —A22- 
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5.2. Differential Equation for [3 



To find a differential equation for (3, into Eq. (5.1) for x put the eigenfunction 

x = 6 exp (iip) 
b = fJ 1 ' 2 



We find then, 



dx 
ds 



db iA 12 , / . / \ 
db + "T" » 6XP W 



using dip/ds = A 12 /b 2 , Eq. (3.15) 



d 



1 ds 



ds \ A\2 ds 



d 



1 db 



+ 



i db 

ds V A12 ds J b 2 ds 
1 iA 12 



db iA\2 
ds b 



A 12 9 



d 



1 db 



ds \A\2 ds 
Putting Eq. (5.4) into Eq. (5.1), one gets 



A 



12 



6 3 



d 



1 db 



exp (iip) 
exp (iip) 
exp (iip) 

A 2 







(5.2) 



(5.3) 



(5.4) 



(5.5) 



ds \ A\2 ds J b 3 ' \ ds \A\2 ) A\2 

Eq. (5.5) is a second order differential equation for b = (3 1 ! 2 . It can be compared to 
the result found when A\2 = 1 and A\\ = 0, 



d 2 b A 



12 



ds 2 6 3 







(5.6) 



5.3. Differential Equation for rj 

r] and x are related by Eq. (4.2) which can be written as 

x = b r), b = (3 1/2 



(5.7) 



In the differential equation for r/ the independent variable is ip or 6 which are related to s 
by 



7 / a ds 
dip = An — 

P 

de = A 12 — R 

u(3 



(5.1 
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We find dx j ds and d(A 1 2 dx / ds) / ds which are then substituted into Eq. (5.1) to get the 
equation for 77, using Eq. (5.5) to eliminate derivatives of b. 



dx db dr/ A\2 



ds 



d_ 

ds 



1 dx 
A\2 ds 



Thus the equation for 17 is 



1 db dr/ A\2 



+ 



ds dip b 
d ( 1 db\ 

ds \A\2 ds J ^ A\2 ds dtp b 2 
d 2 r] A\2 dr/ 1 db 



dtp 2 6 3 
A 12 



+ 



6 3 



d 2 r] A 



+ b 



12 



dtp 2 6 3 
br] \ A21 + 



dtp b 2 ds 
d 

A21 + — 

as 



d ( A u 



An 
A 



12 



+ 



A 2 
A ii 



A 



12 



ds \A 



L12 



+ 



A 
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A 



12 



d 2 r] 
dtp 2 



+ r] = 



d 2 V , 2 

W 2+uri 



The differential equation for r] is unchanged. 



(5.9) 



(5.10) 



5.4. Differential Equations for f3, a, 7 

The differential equations for fl, a, 7, 7 = (1 + a 2 )/ fi can be found starting from Eq. 
(4.13) and (4.14) which we write as 

M = I cos /i + sin /i (5.11a) 
J =(- 7 -a) (5 - Ub) 



We note that 



J(s) =M(s,so)J(so)M(sq,s) 



^-M(s } s ) = A(s)M(s,s ) 
-^M(sq,s) = -M(s ,s)A(s) 



(5.12) 



The last equation follows from M(s, so)M(sq } s) = I. Thus we find that 

d T 

^- = AJ-JA (5.13) 

ds 
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Replacing J using Eq. (5.11b) in Eq. (5.13) we find 

dd 

-£ = (A n - A 22 ) p - 2A 12 a 

n cv 

— = -A 21 (3 - Ani (5-14) 
as 

d^i 

-/ = -2A 21 a- (A n - A 22 )- t 



6. Pertrubation Theory Using the Differential Equation for rj 

The equation for r/ } Eq. (5.10) is often used as a starting point in finding the effects 
of a perturbing field. The particle coordinates are measured relative to a reference orbit 
which is the particle motion in a known magnetic field with components B{. The exact 
equations of motion can then be written as 

^r = J2 Ai i x i + f i z = 1 ' 4 ' -? = 1 ' 4 ( 6 - la ) 

j 

where the fi includes all the terms not included in ^AijXj. These include terms due to 
fields not included in the reference field B{ } which may be referred as AB{ } and nonlinear 
terms due to the terms in the exact equations of motion that do not depend on B{. 

One can see from the exact equations of motion, Eqs. (2.1) that the contributions to 
fi which depend explicitly on ABi, when AB S = 0, are given by 

f2 = ^~(l + x/p) AB y 

B \ (6-lb) 

h = ——( l + x/p) AB X 
Bp 

Repeating the above derivation of Eq. (5.10) for r/ } including the fi terms, one finds the r/ 
equation for the x-motion 

d 2 r] 2 _ vlpl 12 j 
dQ2 + W ~ Au U 



d6 = A 12 - ds 



$x 
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A similar equation can be found for the y motion, 



, R 2 + W = — A — fy 

Jy - U + 1 — 73 + ~r 



A34 (is V^34 

For the case of a gradient perturbation 

A£„ = -Gz (6.4) 

where x = on the unperturbed closed orbit, one can use Eq. (6.2) to find the change in 
u X} Au x . One finds 

If Q 

AUx = 4^ J ds(3x ^ (6 ' 5) 
This well known result for Au x is not changed by using the exact linear equations. 
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